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Abstract. The propagation velocity of the first gas ring in collisional ring galaxies, i.e. the velocity at which the 
maximum in the radial gas density profile propagates radially in the galactic disk, is usually inferred from the 
radial expansion velocity of gas in the first ring. Our numerical hydrodynamics modeling of ring galaxy formation 
however shows that the maximum radial expansion velocity of gas in the first ring (v gaB ) is invariably below the 
propagation velocity of the first gas ring itself (i> r m g ). Modeling of the Cartwheel galaxy indicates that the outer 
ring is currently propagating at « r i ng ~ 100 km s _1 , while the maximum radial expansion velocity of gas in the 
outer ring is currently v gas ~ 65 km s _1 . The latter value is in marginal agreement with the measurements of 
Higdon (1996) based on HI kinematics. Modeling of the radial B — V/V — K color gradients of the Cartwheel 
ring galaxy also indicates that the outer ring is propagating at i> r in g > 90 km s" 1 for the adopted distance to the 
galaxy of 140 Mpc. On the other hand, the azimuthally averaged Ha surface brightness profile of the Cartwheel's 
outer ring does not peak exterior to those in K- and B-bands, contrary to what would be expected for such a 
high propagation velocity. We show that a combined effect of 41° inclination, finite thickness, and warping of the 
Cartwheel's disk might be responsible for the lack of angular difference in the peak positions. Indeed, the radial 
Ha surface brightness profiles obtained along the Cartwheel's major axis, where effects of inclination and finite 
thickness are minimized, do peak exterior to those at K- and 5-bands. The angular difference in peak positions 
implies « r i ng = 110 km s , which is in agreement with the model predictions. We briefly discuss the utility of 
radio continuum emission and spectral line equivalent widths for determining the propagation velocity of gas rings 
in collisional ring galaxies. 
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1. Introduction 

Collisional ring galaxies are thought to originate in galac- 
tic collisions. A near-central axial passage of a companion 
galaxy through the disk of a target galaxy initiates a se- 
ries of expanding ring density waves both in gas and stellar 
components of the target galactic disk. An expanding ring 
gas density wave (hereafter, gas ring) triggers high rates 
of massive star formation along its perimeter, if the gas 
density exceeds a threshold (Appleton & Struck-Marcell 
1987). The propagation velocity of gas ring, i.e. the ve- 
locity at which the maximum in radial gas density pro- 
file propagates radially in the galactic disk, is usually in- 
ferred from the radial expansion velocity of gas in the ring. 
Indeed, kinematical models predict that a determination 
of the radial expansion velocity of an individual star or 
HII region in the first ring gives a good indication of the 
ring propagation velocity (Toomre 1978). This is usually 
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not true for higher order rings. However, rings are rather 
wave phenomena than material features, which implies a 
possible difference between the ring propagation velocity 
and the radial expansion velocity of particles forming the 
ring. Indeed, recent N-body modeling of ring galaxy for- 
mation by Athanassoula et al. Ijl997|) have shown that 
the propagation velocity of the first stellar ring is usually 
higher than the mean velocity of stellar particles in the 
ring. Hence, the propagation velocity of the first gas ring 
may also not coincide with the radial expansion velocity 
of gas in the ring. If so, this would have important im- 
plications for the age estimates of collisional ring galaxies 
(hereafter, ring galaxies), since in many well-known sys- 
tems, most notably in the Cartwheel, the rings are essen- 
tially a gas dynamical phenomenon. 

In this paper we attempt to estimate the propagation 
velocity of gas ring from the optical and near-infrared pho- 
tometry, as well as non-thermal radio emission and spec- 
tral line equivalent widths of ring galaxies. If an expand- 



2 



E. I. Vorobyov and D. Bizyaev: Propagation velocities of gas rings 



ing gas ring indeed triggers high rates of star formation 
as predicted by Appleton & Struck-Marcell l|1987fl . then 
its propagation velocity might be inferred from the prop- 
erties of young stellar component formed in the ring. In 
Sect. |21 we explore the utility of the radial optical and 
near-infrared color gradients measured in the Cartwheel's 
disk (Marcum et al. 1 992) for estimating the propagation 
velocity of the Cartwheel's outer ring. In Sect.|3|we de- 
termine the relationship between the propagation velocity 
of the first gas ring and the maximum radial expansion 
velocity of gas in the ring. The Cartwheel's radial sur- 
face brightness profiles in ifs-band, S-band and Ha are 
compared with those predicted by a density wave model 
in Sect. In Sect. [3] we discuss the utility of radio con- 
tinuum emission and spectral line equivalent widths for 
determining the propagation velocity of gas rings in ring 
galaxies. Our main results are summarized in Sect. El 

2. Radial optical and near-infrared color gradients 
in the Cartwheel's disk 

Theory of ring galaxy formation predicts that in gas- 
rich galactic disks an expanding ring gas density wave 
might trigger a co-existing expanding wave of star for- 
mation (Appleton & Struck-Marcell [1087)- As such a 
density-wave-induced ring wave of star formation prop- 
agates radially from the nucleus, it leaves behind evolved 
stellar populations, with the youngest stars located at 
the current position of the wave. This might result in 
a radial color distribution of stars in the galactic disk, 
with the inner regions being redder than the outer parts 
of the disk. Indeed, observations of the Cartwheel ring 
galaxy (Marcum et al. I1992fl and some other ring galax- 
ies (Appleton & Marston 199TJ) revealed the presence of 
the optical and near-infrared radial color gradients in the 
galactic disks. 

Recent attempts by Korchagin et al. H2()01|l and 
Vorobyov & Bizyaev ({"20011 hereafter VB) to model the 
Cartwheel's radial B — V/V — K color gradients using the 
above scenario yielded a clear success. One of the results of 
Korchagin et al. (2001) was to show that the model radial 
B — V/V — K color gradients are sensitive not only to the 
metallicity of star- forming gas, but also to the propaga- 
tion velocity of gas ring. Indeed, faster rings are expected 
to produce younger stellar populations provided that the 
rings have propagated to the same diameter in the galac- 
tic disk. This motivated us to investigate the utility of the 
Cartwheel's radial color gradients for estimating the ring 
propagation velocity. 

In our previous paper (VB) we postulated the propa- 
gation velocity of the Cartwheel's outer ring as an average 
of all the available measurements (55 km s _1 ) and varied 
the metallicity of star-forming gas in the Cartwheel's disk 
in order to find a better agreement between the model 
and observed radial color gradients. In this paper we per- 
form a more comprehensive modeling of the Cartwheel's 
radial B — V/V— K color gradients and Qbvk color indices 
by varying the propagation velocity as well. Our intuitive 



expectations are that the Cartwheel's radial color gradi- 
ents cannot be reproduced equally well for any value of 
the outer ring propagation velocity. If so, we would be 
able to estimate the outer ring propagation velocity from 
the available optical and near-infrared photometry of the 
Cartwheel ring galaxy. The equations and methods used to 
model the colors and color indices are explained in detail in 
VB and references therein. The population synthesis code 
used in this work is explained in detail in Mayya (1995 
11997(1 . Using stellar evolutionary (Schaller et al. H992J) and 
atmospheric (Kurucz 1992) models, this code synthesizes 
a number of observable quantities in the optical and near- 
infrared passbands, which are suitable for comparison with 
the observed properties of giant star-forming complexes. 
Since the Cartwheel is a recent phenomenon, this code is 
especially suitable for our purposes. The results of this 
code are compared with those of other existing codes by 
Chariot l(T59l)jl . 

A major difference of our modeling from that of VB 
is in the adopted distance to the galaxy, with the value 
of a Hubble constant being Ho — 65 km s _1 Mpc -1 in- 
stead of 100 km s _1 Mpc -1 . The need for this change 
comes from the realization that the results of our mod- 
eling depend on the adopted distance to the galaxy. For 
example, choice of Ho — 50 km s _1 Mpc" 1 instead of 
Hq = 100 km s _1 Mpc -1 would double the linear diam- 
eter of the Cartwheel's disk. The time it takes for the 
density wave of constant propagation velocity to cross 
the Cartwheel's disk would grow accordingly, which will 
inevitably influence the colors of density-wave-born stel- 
lar populations in each radial annulus of the Cartwheel's 
disk. There are distinct signs of convergence between dif- 
ferent estimates of a Hubble constant to the value of 
Ho = 65 ± 10 km s" 1 Mpc" 1 (Binney & Merrifield ITMSjl . 
Therefore, in this work we adopt Ho — 65 km s _1 Mpc -1 , 
which assumes 140 Mpc for the distance to the Cartwheel 
galaxy. 

We consider several propagation velocities of the outer 
ring (fring) starting from 25 km s _1 and ending with 
120 km s _1 . We focus on three values, namely 25 km s~ x , 
55 km s , and 90 km s _1 as implied by the measurements 
of gas expansion velocity by Amram et al. (1998), Higdon 
( 1996), and Fosbury & Hawarden f 19771- For each value of 
propagation velocity we explore a wide range of metallici- 
ties from z — zq/2Q to z = 2zq in order to achieve a better 
agreement between the model and observed photometry. 
The properties of the pre-collision stellar disk are chosen 
typical for the late- type Freeman disks, with the V-band 
central surface brightness [Jy = 21.0 mag arcsec~ 2 and 
the scale length Rq = 3 kpc respectively. In this paper 
we do not model the colors of the Cartwheel's inner ring 
and the nucleus, as it requires additional assumptions on 
the star formation history in the Cartwheel (see VB for 
details) . 

The Cartwheel's and the model colors/indices obtained 
for the outer ring propagation velocity of 25 km s _1 and 
90 km s _1 are shown in Fig.f^ Open triangles in Figs.^, 
and represent the model radial B — V/V — K color 
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gradients and Qbvk combined color indices obtained for 
fring = 25 km s -1 , while the open squares indicate the 
corresponding values obtained for u r i ng = 90 km s _1 . For 
the latter velocity, the model colors of three outermost an- 
nuli are indicated in Fig.^ by Arabic numbers. The filled 
triangles with error bars represent the radial B — V/V — K 
color gradients and Qbvk combined color indices of the 
Cartwheel ring galaxy. The Cartwheel's colors were ob- 
tained by Marcum et al. (1992J for nine annuli and the 
nucleus (see their Plate I) and numbered in Fig. QJl ac- 
cording to radius, beginning with the inner ring (annulus 
I) and ending with the outer ring (annuli VIE- IX). 

The fast outer ring can nicely reproduce the 
Cartwheel's Qbvk indices, if metallicity gradient of star- 
forming gas ranges from z = z Q in the inner parts (annuli 



II-IV) to z 



)/3.7 in the outer parts (annuli V-IX). 



The Qbvk indices minimize the uncertainties in the ob- 
served B — V/V — K colors introduced by dust extinction. 
Hence, a good correspondence between the model and the 
Cartwheel's Qbvk indices indicates that the difference in 
each annulus between the model B — V/V — K colors 
shown in Fig. ^ by the open squares and the Cartwheel's 
B — V/V — K colors might be due to dust extinction in 
the Cartwheel's disk. The value of color excess E(B — V) 
implied by this difference can be used to estimate the 
surface densities of gas in each annulus using standard 
gas-to-dust ratio (Bohlin et al. H978[) . The estimated sur- 
face densities of gas are below the detection upper lim- 
its for each annulus, which indicates that internal extinc- 
tion is indeed responsible for the difference between the 
Cartwheel's B — V/V — K colors and the model colors 
shown in Fig. by the open squares. 

On the other hand, the slow outer ring cannot re- 
produce the Cartwheel's Qbvk indices. Disagreement be- 
tween the model and Cartwheel's indices is less severe, 
if metallicity gradient of star-forming gas ranges from 
z = zq in the inner parts (annuli II-IV) to z = zq/3.7 
in the outer parts (annuli V-IX). Nevertheless, the model 
Qbvk indices shown by the open triangles in Fig. ^> lie 
beyond the observed limits, which indicates that the dif- 
ference in each annulus between the model B — V/V — K 
colors shown in Fig. by the open triangles and the 
Cartwheel's colors cannot be attributed to dust extinction 
in the Cartwheel's disk. 

We considered other propagation velocities of the outer 
ring and found that the Cartwheel's B — V/V — K colors 
and Qbvk color indices could be acceptably reproduced 
for 55 km s _1 < u r j ng < 120 km s _1 . The best-fit value 
is Wring — 90 km s _1 . Stellar populations born in each 
annulus by a slower wave are older than those born by a 
faster wave. Strong dependence of B — V/V — K colors on 
the age of young stellar populations allows us to constrain 
the propagation velocity of the Cartwheel's outer ring. 

Three assumptions can possibly influence our color 
modeling results: the assumed properties of the pre- 
collision stellar disk, the adopted value of a Hubble con- 
stant, and the population synthesis model used in color 
gradient modeling. 




4 5 6 7 

Photometric annulus 

Fig. 1. The radial profiles of a) B — V/V — K color gradi- 
ents and b) the Qbvk combined color indices. The open 
squares and the open triangles show the model profiles ob- 
tained for the outer ring propagation velocity of 90 km s _1 
and 25 km s _1 respectively. The filled triangles with error 
bars indicate a) the measured B — V/V — K radial color 
gradients in the Cartwheel's disk and b) the Qbvk color 
indices calculated from the measured B — V/V — K colors 
of the Cartwheel's disk. The insert shows the results of 
color modeling with Starburst99 (Leitherer et al. [1999). 



a) The pre- collision stellar populations. To understand 
how the model B — V/V — K colors depend on the prop- 
erties of the pre-collision disk, we gradually reduced (Jy 
down to 23.5 mag and increased i?o up to 5 kpc imitating 
transition from the late-type galactic disk to the low sur- 
face brightness (LSB) disk. Again, for v r j ng < 50 km s^ 1 
the model colors and indices totally disagree with those of 
the Cartwheel galaxy. The acceptable agreement between 
the model and the Cartwheel's radial color gradients and 
color indices is found for higher propagation velocities of 
90 km s _1 , though the model Qbvk indices tend to lie sys- 
tematically above the observed values. This again favors 
to the pre-collision Cartwheel being a normal late-type 
galaxy rather than a LSB galaxy (VB). 

b) Distance to the Cartwheel. We recomputed 
the model colors for a Hubble constant of Hq = 
100 km s _1 Mpc -1 , which assumed 90 Mpc for the dis- 
tance to the Cartwheel galaxy. The adopted distance to 
the Cartwheel affects the model radial color gradients. 
A general tendency is that the acceptable agreement be- 
tween the model and the Cartwheel's B — V/V — K colors 
and Qbvk indices can now be achieved for lower values of 
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the ring propagation velocity, more specifically, for v T i ng > 
40 km s^ 1 . The best-fit value is now i> r i ng = 60 km s _1 . For 
the ring propagation velocities higher than 90 km s _1 , the 
correspondence between the model and the Cartwheel's 
colors/indices of the four outermost annuli worsens. 

c) Population synthesis model. We tested our results 
of color gradient modeling with another population syn- 
thesis code, namely Starburst99 (Leitherer et al. |I999). 
This code was developed for synthesizing observable quan- 
tities in starburst galaxies. Since the Cartwheel harbors a 
moderate burst of star formation (SFR w 67 M Q yr _1 ), 
Starburst99 is best suited for comparison with the code 
assumed in this work ( Mayya 1912J). The results of color 
modeling with Starburst99 are plotted in the insert of 
Fig.^, for both the fast outer ring (u r i ng = 90 km s _1 , the 
open squares) and the slow outer ring (t> r ing = 25 km s , 
the open triangles). It is clearly evident that the fast 
outer ring can well reproduce the Cartwheel's colors. On 
the contrary, the slow outer ring cannot reproduce the 
Cartwheel's optical and near-infrared colors. Starburst99 
uses a new atmospheric model of Lejeune et al. CE23 an d 
the latest Geneva evolutionary models, which are different 
from those used by Mayya in his population synthesis code 
(Mayya 1995 ). Nevertheless, the results of color modeling 
with both codes are very similar. Hence, our main conclu- 
sions are not affected by the population synthesis model 
assumed in this work. 

Direct estimate of the Cartwheel's outer ring prop- 
agation velocity from the color gradient modeling per- 
formed in this section is problematic since our results 
are sensitive to the assumed value of a Hubble constant. 
Nevertheless, our modeling indicates that the propaga- 
tion velocity of the Cartwheel's outer ring cannot be lower 
than 40 km s _1 for any reasonable value of a Hubble con- 
stant. Determining the ring propagation velocity may be 
in fact even more difficult, because it may not coincide 
with the expansion velocity of gas in the ring. Indeed, 
rings are rather wave phenomena than material features 
(Athanassoula et al. H997fl and the difference between the 
ring propagation velocity and the radial expansion veloc- 
ity of particles forming the ring is not totally unexpected. 
In the next section two dimensional numerical simulations 
of a Cartwheel-like ring galaxy are performed. We restrict 
ourselves to the planar geometry in order to find the rela- 
tionship between the maximum radial expansion velocity 
of gas in the ring and the propagation velocity of the gas 
ring itself. 

3. Numerical simulations of a Cartwheel-like ring 
galaxy 

3.1. Initial conditions 

An extensive study of the Cartwheel galaxy by Higdon 
(1995, 1996) provided important observational constrains 
on the parameters of our numerical model. The indicative 
mass of the Cartwheel, which generally represents the to- 
tal mass to within a factor of about 2, is 5.3 x 10 11 M m 



for the adopted H — 65 km s _1 Mpc -1 (Higdon 
An estimate of the luminous stellar mass M st in the 
Cartwheel is derived by multiplying its blue luminos- 
ity L B = 3.7 x 10 10 Lq (Higdon ITMty by a late spi- 
ral ratio M/L B =3 (Faber & Gallagher EHZOJ- This re- 
sults in = 1.1 x 10 11 M Q . The total HI mass of 
the Cartwheel is 2.2 x 10 10 M fffigdon [T95SI corrected 
for Ho = 65 km s _1 Mpc -1 ), with the total gas mass 
(M m + M He ) amounting to 3 x 10 10 M Q . The Cartwheel 
galaxy thus appears to be halo-dominated. This assump- 
tion was also used by Appleton & Higdon (1993) in their 
hydrodynamic models of the Cartwheel galaxy. 

We consider a Cartwheel-like galaxy consisting of a 
sclf-gravitating thin gaseous disk of M& = 3 x 10 10 M and 
a rigid halo of M h = 5 x 10 11 M . Marcum et al. I|1992|l 
found no evidence for old pre-collision stellar populations 
in the outer parts of the Cartwheel. Numerical model- 
ing of the radial B — V/V — K color gradients confirmed 
that the pre-collision stellar populations dominate only 
in the inner parts of the Cartwheel (VB). Moreover, our 
numerical simulations indicate that stellar density waves 
propagate only 2-3 disk scale lengths and become signifi- 
cantly less defined as compared to gas density waves, the 
result also found in numerical simulations by Hernquist 
& Weil (1993). Hence, we expect that the stellar density 
wave would have little effect on the dynamics of gas at the 
radii of the Cartwheel's outer ring. 

The radial distribution of gas in the pre-collision 
Cartwheel is unknown. Our modeling of radial color gra- 
dients measured in the Cartwheel's disk (VB) suggested 
that this distribution most probably showed a mild expo- 
nential decline with a surface density contrast between the 
nucleus and the position of the outer ring being a factor of 
3.5. More importantly, the lack of old stellar component at 
the outer parts of the Cartwheel implies that the surface 
density of gas in the pre-collision Cartwheel was below the 
critical surface density for gravitational instability every- 
where except probably the innermost radii. The critical 
surface density is defined as 
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where v s is the speed of sound and k is the epicyclic fre- 
quency (Martin & Kennicutt l2TM)l . 

Initially, the rotating gaseous disk is balanced by self- 
gravity of the disk, the radial pressure gradient, and grav- 
ity of the spherically symmetric rigid halo 
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dr 



(2) 



v^ = J_dP 

r £ g dr 

where v c i T is the rotational velocity, P is the gas pressure, 
$h is the halo potential, and $d is the potential of gaseous 
disk defined in the polar coordinates as 
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We adopt an isothermal equation of state (T = 
10 4 K°), which is an usual approximation used for mod- 
eling the interstellar medium in disk galaxies (Mihos & 
Hernquist I1994[) . Struck (|1997[) has demonstrated that 
the ring structure and dynamics in the isothermal and 
non-isothermal simulations are similar. We also tried an 
adiabatic equation of state and found little difference in 
the first ring structure and dynamics, though the second 
ring structure and dynamics showed substantial differ- 
ences with respect to the isothermal case. 

A slightly modified form of rotation curve in a soften- 
ing point-mass potential is adopted 



c(0 



v r 



.0.9 



(r 2 + i?2)0.5 ' 



(4) 



where Rq is the effective halo radius. The power index 0.9 
allows us to model a Keplerian decline in rotation curves 
at larger radii. Parameters of the rotation curve (eq. 0J) 
are chosen so that the halo mass 



M h = 



G dr 



-Rou 



(5) 



is equal to 5 x 10 11 M . As the rotational curve is fixed by 
eq. and the disk potential $d is determined by numer- 
ically integrating eq. J3J), the gradient of halo potential 
$h can be calculated using eq. (J2J. Parameters used in 
numerical simulations are listed in Tabled 

Table 1. Model parameters of a Cartwheel-like galaxy 



Parameter 


Value 


Inner disk radius (Ri n ) 


0.2 kpc 


Outer disk radius (Rout) 


34 kpc 


Disk scale length (ro) 


20 kpc 


Central surface density of gas (E g (0)) 


23.5 M Q pc~ 2 


Gas temperature (T) 


10 4 K° 


Disk mass (Ma) 


3 x 10 10 M @ 


Halo mass (Mh) 


5 x 10 11 M© 


Mass of a companion (M c ) 


1.5 x 10 11 M @ 


Companion's softening parameter (e) 


5 kpc 



3.2. Computational techniques 

We use the ZEUS-2D numerical hydrodynamics code in- 
corporating a self-gravitating thin gaseous disk and a 
fixed spherically symmetric halo potential. A usual set of 
isothermal hydrodynamics equations in the polar coordi- 
nates (r, cj>) is solved using the method of finite-differences 
with a time-explicit, operator split solution procedure ex- 
plained in detail in Stone & Norman (1992). The inter- 
action is treated by computing the gravitational force ex- 
erted on the disk by the softened point-mass potential of 
a companion galaxy. Equations of motion of the compan- 
ion in the combined gravitational potential of the halo 
and the gaseous disk are numerically solved using the 
Runge-Kutta scheme. The gravitational potential of the 



thin disk (eq. |3J) is evaluated numerically using the fast 
Fourier transform (Binney & Tremaine 1987). The code 
was verified against a standard set of test problems de- 
scribed in Stone & Norman $T$92). 

3.3. Results of modeling. 

We attempt to find the relationship between the propaga- 
tion velocity of ring gas density wave and the maximum 
radial expansion velocity of gas in the wave generated by 
a near central passage of a companion galaxy through the 
disk of a Cartwheel-like galaxy. We have chosen the model 
parameters typical for the Cartwheel galaxy in order to 
compare the results of numerical simulations with those 
of Sect. [3 By varying the initial rotation curve, the mass 
of a companion, and the impact point we try to repro- 
duce the morphology, kinematics, and radial distribution 
of gas in the Cartwheel's disk. Below we briefly describe 
the parameter space investigated in our simulations. 

Initial rotation curve. We consider three initial rota- 
tion curves shown in Fig. [21 and find that the Cartwheel's 
HI morphology and kinematics are best reproduced for 
A2 curve. Al curve produce too small a spacing between 
the inner and outer rings as compared to that of the 
Cartwheel, while A3 curve does not develop the inner ring 
at all. 

Mass of a companion. We consider a wide range of 
companion galaxy masses from 4x 10 10 M Q to 3x 10 11 Mq, 
which is roughly 7.5% — 55% of the Cartwheel's total mass. 
The indicative mass of a companion G3, most massive 
in the Cartwheel's group (Hiedon H"§96). is A/q3 = 4 x 
10 10 Mq. This mass is not enough to reproduce the radial 
surface density profiles of HI in the Cartwheel. Such a low- 
mass companion gives a weak first response, with less than 
a 20% enhancement in gas densities over the unperturbed 
values. At the time when the first ring reaches 24 kpc 
radius (the current location of the Cartwheel's outer ring) , 
the actual mass of gas locked in the ring is only 25% of 
the total gas mass, while in the Cartwheel the outer ring 
contains about 85% of total HI (Higdon ITMfJ|l . The low 
mass companion does not reproduce the observed spacing 
between the inner and the outer Cartwheel's rings, with 
the radius of the model second ring being almost twice 
bigger than that of the Cartwheel's inner ring. 

We find that companions with masses no less than 25% 
of the Cartwheel's total mass can acceptably reproduce 
the HI distribution in the Cartwheel's disk. In our simu- 
lations the companion's mass is 1.5 x 10 11 Mq, which is 
roughly 30% of the total mass of the Cartwheel. A 30% 
companion generates stronger response, with more than a 
100% enhancement in gas surface densities over the unper- 
turbed values. There are no companions to the Cartwheel 
large enough to fit the companion's mass we determine. 
The companion's optical and HI line widths give mass es- 
timates that are only ~ 6% of the Cartwheel's (Davies 
& Morton IT982I Higdon ESEJl • However, determining ac- 
curate galaxy masses is difficult enough in strongly per- 
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Fig. 2. Rotation curves of pre-collision gaseous disk gener- 
ated using eq. 0] Parameters determining rotation curves 
are: v Q = 375 and R = 4 (marked "Al"), v = 382 
and R = 10 (marked "A2"), v = 420 and R = 18 
(marked "A3"). The Cartwheel's outer ring is at 24 kpc 
(Ho = 65 km s" 1 Mpc^ 1 ). 

turbed galaxies like the Cartwheel's companions. Passing 
through the Cartwheel's disk causes the companion to lose 
mass due to tidal and ram pressure stripping. Indeed, 
Higdon ( 1996) reported HI plumes stretching from the 
Cartwheel toward G3 companion. Ring galaxy simulations 
by Tsuchiya et al. 1)199 8?) found that massive companions 
(25% of ring galaxy's mass) can loose around 35% of their 
halos, or nearly a third of their total mass, due to the in- 
teraction. The gas disks in massive companions were also 
greatly truncated (R gas < 50% of original radius). If much 
of the companion's mass resided in a halo, one would in 
fact significantly underestimate its total mass. The appar- 
ent low masses of the companion galaxies are a problem, 
but it is likely to be a consequence of the ring galaxy for- 
mation process. More work is clearly needed in estimating 
masses of ring galaxy companions to see if this is a general 
result. 

Impact points, impact angles, and relative velocities. 
In our simulations the companion placed initially at a 
distance of 35 kpc above the Cartwheel's plane, which 
was roughly 150% of the Cartwheel's radius, was given 
an initial relative velocity of 80 km s _1 toward the 
Cartwheel. Varying the initial distances from the axis 
of the Cartwheel, we obtained a wide range of im- 
pact points starting from a nearly bull's-eye collision 
(-Rimpact=0.1 kpc) and ending with an essentially off- 
center collision (-Ri m pact=10 kpc). The former produces 
symmetric rings with negligible azimuthal variations in 
gas surface densities around the ring, while the latter re- 
sults in crescent or spiral-like density waves. Neither of 
these ultimate cases resembles the actual HI distribution 
in the Cartwheel's outer ring. We find that the Cartwheel's 



See figure3.gif 

Fig. 3. The model velocity field of gas superimposed on 
the contour plot of gas surface density £ g obtained for 
A2 curve and model parameters listed in Table H The 
solid line outlines the regions where £ g exceeds the critical 
surface density for gravitational instability defined by eq. 
n The scale bar is in M Q pc -2 . 

HI kinematics and morphology are best reproduced for 
a slightly off-center i?i m pact=2.5 kpc collision. The rela- 
tive velocity of the companion at the time of impact is 
480 km s" 1 . We also considered other distances above the 
Cartwheel's plane and other relative velocities as the ini- 
tial parameters of the companion. For the initial distance 
of 70 kpc and initial relative velocity of 80 km s _1 toward 
the Cartwheel, the relative velocity at the time of impact 
was Relative (impact) = 540 km s _1 . For the initial distance 
of 35 kpc and initial relative velocity of 380 km s _1 toward 
the Cartwheel, u re iative (impact) was 610 km s _1 . Hence, 
^relative (impact) is largely determined by the parameters 
of the Cartwheel. Variations of 15% in the relative velocity 
affect insignificantly the gas dynamics after the collision. 
Angles of impact from the axis of the Cartwheel are 5-15° 
in all considered cases. We have not made special efforts 
to model highly inclined collisions. 

Figure shows the model velocity field of gas super- 
imposed on the contour plot of gas surface density £ g 
obtained for A2 curve and model parameters listed in 
Table ^ The contrast between £ g in the outer ring and 
the inter-ring region is in agreement with observations. 
S g between the rings in Fig. [3] never exceed 7 M pc~ 2 
and falls below 5 M Q pc~ 2 behind the densest regions 
in the outer ring. Maximum and average £ g around the 
outer ring are 18.5 and 15 M Q pc -2 , however the ob- 
served high density condensations up to 60 M Q pc -2 are 
not found in simulations. E g is below 7 Mq pc -2 beyond 
the outer ring. High density clumps are found in the in- 
ner ring and central regions, with S g « 60-100 M pc~ 2 
and even higher. Development of such high density clumps 
is interesting, since high densities and low temperatures 
are thought to favor HI — > Hi transitions in giant gas 
clouds. Indeed, Horrcllou ct al. (1998) detected a substan- 
tial amount of molecular hydrogen in the Cartwheel's in- 
ner part. For the solar metallicity implied by color gra- 
dient modeling in Sect. [21 Hi mass within the inner ring 
amounts to 2.1 x 10 9 M© according to Horellou et al. 
(1998). In our simulations, the inner ring contains < 10% 
of the gas mass or < 3 x 10 9 M©, which is roughly consis- 
tent with the detected Hi mass. In Sect. I3.ll we neglected 
a possible contribution of molecular hydrogen to the gas 
mass of the Cartwheel as it was only a small fraction to 
the total mass of the Cartwheel galaxy (< 1%) and would 
not affect noticeably the companion's dynamics. 

While the inner ring in Fig. EH has on average higher 
surface densities of gas than the outer ring, it contains 
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Fig. 4. The propagation velocities of the outer ring (the 
filled squares) and the maximum expansion velocities of 
gas in the ring (the open circles) traced along four different 
polar angles as indicated in each frame. The propagation 
velocities of the outer ring are systematically higher than 
the maximum expansion velocities of gas in the ring. 



less than 10% of the total disk mass. On the contrary, the 
outer ring in Fig. [21 comprises more than 50% of the to- 
tal disk mass. No obvious spokes were obtained for either 
of three initial rotation curves. Indeed, HI observations 
of Higdon (1996) have shown that the Cartwheel's spokes 
are HI poor. The solid line in Fig. [2] outlines the regions 
where £ g exceeds the critical surface density for gravi- 
tational instability defined by eq. l{T|l. The regions where 
star formation is supposed to occur correlate with the re- 
gions of maximum £ g in the outer ring. This is not the 
case for the Cartwheel, where a pronounced anticorrela- 
tion between Ehi and Ha surface brightness was found 
(Higdon 1996). Among other reasons such as a blow-out 
in giant supershells formed by multiple supernova explo- 
sions, this might be explained as the consumption of the 
gas supply by the starburst. Inclusion of star formation 
processes in the code is required to check this hypothesis. 

The velocity field in Fig. [21 shows all basic features of 
HI kinematics found in the Cartwheel: (1) a general in- 
crease in expansion velocity of gas for R > 16 kpc, (2) 
inflow of gas towards the inner ring for R < 16 kpc. The 
flow of gas inside the inner ring is complicated. However, 
a general expansion of gas is clearly seen in this innermost 
region. We determine the relationship between the propa- 
gation velocity of the outer ring w r i n g (i.e. the propagation 
velocity of gas density wave) and maximum radial expan- 
sion velocity of gas w gas in the outer ring at different ring 
expansion phases. The regions of maximum density en- 
hancements in the outer ring are used to trace the outer 
ring position during the runs. We start from very early 
times when the outer ring was only rj 6 kpc in radius and 
stop at the outer ring current position of ~ 24 kpc. Due 



to obvious asymmetric structure of £ g in Fig. [31 these ve- 
locities were traced along four different polar angles of 0°, 
90°, 180°, and 270°. The results are shown in Fig.H The 
propagation velocities of the outer ring (the filled squares) 
are systematically higher than the maximum radial expan- 
sion velocities of gas in the ring (the open circles). Both 
velocities generally decline as the outer ring propagates 
outwards. At the time when the outer ring reaches its 
current location at the radius of 24 kpc, w r ing and i> gas 
slow down on average to 102 km s _1 and 65 km s _1 re- 
spectively. However, sinusoidal azimuthal variations in the 
values of both velocities are evident in Fig. 01 The value 
of w r i n g = 102 km s _1 is consistent with the estimates of 
Sect. El where best agreement between the the model and 
observed radial B — V/V — K color gradients was found 
for v r i ng > 90 km s _1 . The value of v gas = 65 km s" 1 
is close to the upper estimate of Higdon (1996 ) based on 
HI kinematics in the Cartwheel's outer ring. The value of 
w gas = 13 — 30 ± 10 km s _1 obtained by Amram et al. 
(1998) using Ha kinematics of ionized gas stands apart. 
It is known that dynamics of different interstellar phases 
may considerably differ in disk galaxies. Inclusion of mul- 
tiphase gas dynamics and star formation in the code is 
needed to determine the dynamical properties of ionized 
gas in ring galaxies. 

We performed a series of computations with different 
values of M n , M c , Tq, and different initial rotation 
curves. In Figs, [5M5b we plot the maximum expansion 
velocities of gas in the ring at the radii of 12 (the up- 
per frame), 18 (the intermediate frame), and 22 kpc (the 
lower frame) as a function of the azimuthally averaged 
propagation velocities of the first ring at the same radii. 
Parameters of each run labeled by Arabic numbers in 
Fig. are listed in Table [2 Impact points and impact 
angles vary slightly in the range 1.5-3 kpc and 5-15° re- 
spectively. We find that the inequality w r i ng > w gas holds in 
all cases, though the resulted morphology and kinematics 
of gas may hardly resemble those of the Cartwheel. Both 
velocities are generally higher at the smaller radii and they 
decline at the larger radii. Higher mass companions pro- 
duce a stronger response with higher ring propagation ve- 
locities and maximum radial expansion velocities of gas 
in the ring, thus leaving the relation between them un- 



changed. It appears that the v T 



> v, 



gas 



inequality holds 



for the first gas rings triggered by the passage of a compan- 
ion through the disk of a halo-dominated target galaxy. 
The least-squares fit to a linear function w gas = k w r ; ng 
shown in each panel of Fig. by a dash-dotted line gives 
k = 0.67 for the rings of 12 kpc, k = 0.55 for the rings of 
18 kpc, and k = 0.51 for the rings of 24 kpc. A linear func- 
tion u gas = Wring is plotted for comparison in each panel 
of Fig. by a dashed line. We also considered gas rings in 
ring galaxies with comparable disk and halo masses, for 
which we obtained an average value of k = 0.8 for the 
rings of 12, 18, and 24 kpc. In higher order resonant rings 
(i.e. the second rings) the flow is more complicated and 



> v, 



gas 



inequality is no more valid. The v- 



ring 



> v e 



inequality contradicts the predictions of a simple kinemat- 
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Fig. 5. The maximum expansion velocities of gas in the 
ring as a function of propagation velocities of the first 
ring at the radii of a) 12 kpc, b) 18 kpc, and c) 22 kpc. 
Parameters of nine runs labeled by Arabic numbers in the 
upper frame are listed in Table [21 The least-squares fit to 
a linear function u gas = k v T i ng is shown in each panel by 
a dash-dotted line. The coefficient k is defined in the text. 
A linear function v gas = v r i ng is plotted for comparison in 
each panel by a dashed line. 



ical model of ring galaxy formation, which argues that a 
determination of the expansion velocity of an individual 
star or HII region in the first ring gives a good indication 
of the ring propagation velocity (Toomre ll978[) . 

Higdon ( 1996) reported the narrow HI line widths and 
the absence of strong shear in the Cartwheel's outer ring. 
He concluded that the outer ring appeared to initiate star 
formation activity primarily through concentration of the 
interstellar medium in the ring rather than through high- 
speed disruptive cloud collisions (^relative > 100 km s~ x ) 
as advocated by Olson & Kwan (1990) for starburst ac- 
tivity in interacting/merging galaxies. Apparent lack of 



molecular hydrogen in the outer ring (Horellou ll995[) ar- 
gues against compact dense clouds. Hence, we have cho- 
sen a smooth initial gas distribution for the pre-collision 
Cartwheel rather than a clumpy medium. In the wider 
context of a disparity in ring propagation velocity versus 
gas expansion velocity in ring galaxies, the i> r i ng > v gas 
inequality has to be proved for a clumpy gas distribution. 
In this connection, sticky particle codes (Combes & Gerin 
1985 ) could be very useful. We note here that N-body sim- 
ulations of stellar rings in ring galaxies by Athanassoula et 
al. (1997) have also shown that the expansion velocity of 
the first stellar ring is usually higher than the mean veloc- 
ity of stellar particles constituting the ring (their Fig. 17, 
the upper frame). Hence, a disparity in ring propagation 
velocity versus expansion velocity of gas/stars constitut- 
ing the ring might be a general property of ring galaxies. 



Table 2. Model parameters of nine runs shown in FigEl 
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The w r i ng > fgas inequality for the first gas ring in ring 
galaxies has an important implication. Measurements of 
the radial expansion velocity of gas in the Cartwheel's 
outer ring based on HI and Ha kinematics appear to 
underestimate the true propagation velocity of the outer 
ring. Moreover, the radial expansion velocity of gas itself 
may be underestimated. Indeed, w gas is derived assuming 
that all the gas motions are in the plane of the ring. The 
plane geometry neglects any non-planar component of v ga , s 
that may be present in warped rings. Results of numeri- 
cal experiments performed by Struck (1997) indicate that 
the Cartwheel's disk might be slightly warped. Hence, the 
measurements of gas expansion velocity in the Cartwheel's 
outer ring may underestimate v gas , which in turn under- 
estimates the true propagation velocity of the outer ring. 
The same might be true for other ring galaxies. If so, ring 
galaxies are generally younger than they have been previ- 
ously thought. If rings are younger then we can expect to 
see even fewer of them - they should be relatively rarer. 
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See figure3.gif 

Fig. 6. The model Jf-band (filled squares), 5-band (open 
circles), and Ha (filled triangles) radial surface bright- 
nesses of the Cartwheel ring galaxy obtained for the same 
model parameters as in Fig.Q] The surface brightness pro- 
files peak at/near the position of the outer ring at 24 kpc. 
Difference in the peak positions at if-band and Ha is big- 
ger than at 5-band and Ha. 



4. Radial 5-band, if-band, and Ha surface 
brightness profiles of the Cartwheel galaxy 

Marston & Appleton l|1995[) found that unlike normal 
galaxies the azimuthally averaged radial intensity profiles 
in Ha and JT-band of most ring galaxies peak at the po- 
sition of the ring. Specifically, for the majority of ring 
galaxies the peak in the Ha radial profile lies at the larger 
radii than the peak in the if-band radial profile. Indeed, 
a standard scenario of ring galaxy formation predicts that 
the if-band and Ha peaks in radial intensity profiles of 
density-wave-born stellar populations are spatially sepa- 
rated, with the peak in Ha lying on the leading edge of 
the wave where most star formation takes place, and the 
peak in if-band lying progressively behind the wave as the 
density-wave-born stellar populations evolve (Korchagin 
et aims)- 

Applying the model used to compute radial color gra- 
dients in Sect [3 we found that the linear difference be- 
tween the if-band and Ha peak positions depends not only 
on the wave velocity, but also on the metallicity of star- 
forming gas. Specifically, this difference is proportional to 
the wave velocity. For metallicities of star-forming gas be- 
low 1/5 of the solar there is no peak in iT-band at the 
position of the ring and the model if-band radial sur- 
face brightness profile shows a monotonous decline at all 
radii. The peak appears at higher metallicities and be- 
comes more pronounced with increasing metallicity. In 
Fig. we present the model azimuthally averaged radial 
surface brightness profiles of the Cartwheel galaxy in K- 
band (filled squares), B-band (open circles), and Ha (filled 
triangles) . All model parameters of the gas density wave 
and pre-collision stellar populations are chosen as in Fig.n 
The model surface brightness profiles peak at/near the po- 
sition of the Cartwheel's outer ring at 24 kpc. Analysis of 
Fig. indicates that the difference between the Ha and 
if-band peak positions is expected to be around 1.1 kpc 
(1.65") for the adopted distance of 140 Mpc and the outer 
ring propagation velocity of 90 km s _1 . The Ha — B peak 
difference is 0.7 kpc (1.05"). We note here that the un- 
derlying ring stellar density wave may influence the peak 
position in 7f-band. To exclude this effect, density- wave- 
born stellar populations must dominate the ring emission, 
which appears to be the case for the Cartwheel galaxy 
(Marcum et al. lTM21 VB). 



See figure3.gif 

Fig. 7. The Cartwheel galaxy in Ks-b&nd. The dashed 
ellipses outline elliptical annuli used for determining the 
azimuthally averaged radial surface brightnesses of the 
Cartwheel. The north is shown by the arrow. 

We compare our model predictions with the actual 
angular differences between the peaks in the azimuthally 
averaged 5-band, JTs-band and Ha radial surface bright- 
ness profiles of the Cartwheel's outer ring. Images of the 
Cartwheel together with the calibrating frames were ob- 
tained from the public archives. The near-infrared images 
in ifs-band were provided by ESO VLT archive. Images 
in the HST's WFPC2/F450W (B) band were taken from 
the Hubble Space Telescope archive. Hereafter, we refer 
to this band as the "5-band". Images in Ha and "Ha- 
free continuum" were provided by the CFHT archive. All 
the data were reduced in the standard way. Since we were 
interested only in determining the angular difference be- 
tween the peak positions in different passbands, the flux 
calibration was not applied. All the images were aligned 
using the standard MIDAS routines. The estimated accu- 
racy of alignment is better than 0.03". Finally, an angular- 
resolution of 0.76" was achieved for the Ha image. When 
the B and Ks images are considered separately from the 
Ha image, a better resolution of 0.46" is achieved. At the 
adopted distance of 140 Mpc, this corresponds to a spatial 
resolution of 500 pc for the Ha image and 300 pc for the 
images in Ks- and 5-bands. 

Figure [7| shows the image of the Cartwheel galaxy in 
Ks-band. The north is down and the east is to the right. 
As is seen in Fig. the geometrical center of the outer 
ring does not coincide with the nucleus. To take this pecu- 
liar morphology into account, a series of elliptical annuli 
were constructed, with the smallest ellipse centered on the 
nucleus and the centers of larger ellipses slightly shifted so 
as the largest ellipse is centered on the geometrical center 
of the outer ring. These elliptical annuli were then used 
to determine the azimuthally averaged ifs-band, 5-band, 
and Ha radial surface brightness profiles shown in Fig. |SJ 

Surprisingly, we find no noticeable angular difference 
between the peak positions of the azimuthally averaged 
radial surface brightness profiles in the Cartwheel's outer 
ring, contrary to what is expected from our modeling (see 
Fig. EJ). Absence of such a difference is also evident in 
Fig. 4 of Higdon (1993- This implies that either the outer 
ring velocity is too slow for this difference to be resolved or 
considerable inclination of the Cartwheel acts to smear out 
this difference. Indeed, if the Cartwheel's disk is warped 
("Struck fl997|l . then the elliptical annuli might not be a 
good choice for obtaining the azimuthally averaged radial 
profiles. Another geometry effect may also affect the az- 
imuthally averaged profiles in considerably inclined ring 
galaxies. For an inclined disk of finite thickness, the line- 
of-sight flux comes from density-wave-born stellar popu- 



10 



E. I. Vorobyov and D. Bizyaev: Propagation velocities of gas rings 



u 

CD 
00 

u 



cn 

a 

E 

uf 
m 

CD 

c 



CD 

CD 
U 

a 



20 



22 \ 



24 



26 





K + - 


- - + - 




B - 




V 


Ha 


— A 

- 


\ \ 




- 

r — 


IN. 


i 
f 

/ 


\ - 




/ A 


|V - 










N 


\ 

I 


V 








\ 






l 


1 1 




i 

i 

i 



10 20 30 40 
Radius, arcsec 

Fig. 8. The azimuthally averaged surface brightness pro- 
files of the Cartwheel galaxy. All profiles peak at the po- 
sition of the outer ring at 36". No difference between the 
peak positions is seen. 



lations located at different radii, and hence having differ- 
ent ages and different photometric properties. This line- 
of-sight overlapping of stellar populations of different ages 
may have a significant effect on the azimuthally averages 
surface brightness profiles, if most of the radiation flux 
comes from regions close to the minor axis and the line- 
of-sight thickness of the disk is comparable to the expected 
shift between the peaks in different passbands. 

Rough estimates of this effect can be done by assum- 
ing that the total stars+gas surface density in the ring 
(Stot) is of order 20 M© pc -2 and using a typical value of 
the stellar velocity dispersion of young stellar populations 
c z = 10 km s _1 . For the case of a self-gravitating disk, 
the isothermal scale height of the ring h = c 2 /27rG£tot is 
near 0.2 kpc. The vertical thickness of the outer ring 2h 
is thus 0.4 kpc. The gravitational potential of a dark halo 
makes a disk thinner, but not very significantly (Bahcall 
1984J). At the same time, disk thickness is expected to 
be twice bigger for interacting galaxies (Reshetnikov & 
Combes Il997fl . Results of numerical experiments also in- 
dicate that the Cartwheel's outer ring might be thick and 
slightly warped (Struck IT537|l . 

Our model predicts that the linear difference between 
the Ha and JT-band peaks is of order of a kiloparsec and 
this value is thus comparable to the line-of-sight thickness 
of the Cartwheel outer ring for the adopted inclination of 
41° (Higdon 199(f). Hence, we expect that both the line- 
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Fig. 9. The angular difference between the peaks in radial 
surface brightnesses azimuthally averaged over a 1° sector 
as a function of position angle measured counterclockwise 
from a line extending north from the outer ring center, a) 
Angular difference between the Ha and ifs-band peaks, 

b) angular difference between the Ha and 5-band peaks, 

c) angular difference between the Ks- and i?-band peaks. 
Solid lines indicate the angular peak differences smoothed 
by a wide median filter (30°). 



of-sight overlapping and the warping act to smear out the 
difference between the peak positions of the azimuthally 
averaged radial surface brightness profiles shown in Fig. |SJ 
Fortunately, overlapping should have little effect in the 
regions located along the major axis of the Cartwheel. 

We apply an alternative method for determining the 
angular difference between the peak positions of radial 
profiles in different photometric bands. Instead of com- 
puting the fluxes from each elliptical annulus shown in 
Fig- El i- e - t ne fluxes averaged over the whole azimuth, we 
make a series of radial cuts drawn from the Cartwheel's 
nucleus with an angular separation of 0.2°. A series of 360 
radial surface brightness profiles of the Cartwheel's outer 
ring were then constructed by azimuthally averaging these 
cuts over a 1° sector. Peaks were located by both fitting 
a Gaussian to each profile and determining the "center of 
mass" of each profile. The radial positions of each peak 
were corrected for a 41° inclination of the Cartwheel and 
re-calculated for a zero inclination. 
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Both methods give very similar results. Angular dif- 
ferences in arcsec between the peak positions are shown 
in Fig. El as a function of position angle, which is mea- 
sured counterclockwise from a line extending north from 
the outer ring center. The ifs-band and 5-band peaks are 
located at smaller radii than the peak in Ha for almost all 
position angles. The angular difference between the Ks- 
band and Ha peaks is systematically bigger than between 
the -B-band and Ha peaks (Fig. Et>), which is in full agree- 
ment with our model predictions. 

Let us now consider the effects of finite thickness of 
the Cartwheel's disk in more detail. Let us assume that 
the Cartwheel's disk has a radially non-constant thickness, 
with the inner regions being thicker than the outer regions. 
This assumption is justified, since the inner regions of the 
Cartwheel are expected to be older and hence thicker than 
the outer regions. Internal dust extinction is known to 
be considerable in the Cartwheel's outer ring (Fosbury & 
Hawarden 119771 VB). We assume a limiting case when 
part of the Cartwheel's outer ring that is located below 
the galactic plane is totally unseen. Under these assump- 
tions, the angular differences between the peak positions 
should have maximums along the minor semi-axis that is 
drawn on the galaxy's near side. The minima of angular 
differences are then expected to be along the minor semi- 
axis drawn on the galaxy's far side. We smoothed Fig. 
with a 30° wide median filter in order to minimize the 
noise. The results are shown by the solid lines in each 
panel of Fig. [5J The angular differences between all three 
bands are minimal along the southern minor semi-axis at 
position angles 180°-220°. On the opposite side, along the 
northern minor semi- axis at position angles 20°-40°, the 
angular differences in Fig. [5| show an expected tendency 
to grow and reach a maximum. This places the southern 
part of the Cartwheel on the galaxy's far side. Our con- 
clusion is consistent with the analysis of Higdon (1996), 
who deduced the galaxy's orientation indirectly from the 
form of rotation curve and spokes. There is however an 
unexpected maximum in the Ha — K and Ha — B an- 
gular peak differences at position angles 320°-340°, while 
the B — K angular peak difference is minimal there. This 
may be a manifestation of an Ha warp in the Cartwheel's 
outer ring. 

As it has been noted above, differences between 
the peak positions measured along the major axis are 
not influenced by the finite thickness of a ring galaxy. 
Hence, they might be considered as an indicator of the 
Cartwheel's outer ring propagation velocity. The angu- 
lar peak differences along the Cartwheel's major axis de- 
termined from the smoothed curves in Fig. [S] are Ha — 
Xs=1.97", Ha - 5=1.29", and B - Ks=0.67". They cor- 
respond to the linear peak differences of Ha — Ks= 1.3 kpc, 
Ha — i?=0.85 kpc, and B — Ks— 0.45 kpc for the adopted 
distance of 140 Mpc. These values are even bigger than 
those predicted by our modeling for the Cartwheel's outer 
ring velocity of 90 km s _1 , implying that the outer ring 
is currently propagating at w r i ng « 110 km s . We note 
that these estimates cannot be considered as very accu- 



rate, though they clearly show that the Cartwheel's outer 
ring is indeed expanding, and probably with a consider- 
able velocity as predicted by modeling of color gradients in 
Sect. Eland numerical hydrodynamics modeling in Sect.|3 
We derive v r i ng « 100 km s _1 as an average of estimates 
obtained in Sect. 12141 Since the wave velocity most prob- 
ably decrease with time, we conclude that the Cartwheel 
galaxy is no older than 250 Myr. 

5. Other possibilities for estimating the 

propagation velocity of ring gas density waves 

5.1. Call triplet equivalent widths 

The Call triplet (CaT) lines originate in the atmospheres 
of cool stars with the absorption being the strongest in red 
supergiants. While studying numerically the evolutionary 
behavior of the CaT equivalent widths (EW(CaT)) in star- 
bursts of different metallicities, Mayya (1997) found that 
besides the primary peak around 10 Myr the EW(CaT) 
show a well defined secondary peak around 60 Myr in the 
solar metallicity starbursts. The time separation tcaT = 
50 Myr between the primary EW(CaT) peak around 
10 Myr and the secondary EW(CaT) peak around 60 Myr 
might be in principle used to infer the propagation veloc- 
ity of gas density waves in ring galaxies. 

We illustrate this possibility for the Cartwheel ring 
galaxy. We assume that the Cartwheel's outer ring is cur- 
rently propagating at v T i ng =90 km s . For the time be- 
ing, we assume that the Cartwheel's gaseous disk has 
the solar metallicity. We also neglect the possible input 
to EW(CaT) from the pre-collision stellar populations, 
which is indeed minimal at the Cartwheel's outer parts 
(VB). The filled triangles in Fig. 1101 show the model ra- 
dial profile of EW(CaT) after the outer ring has propa- 
gated 24 kpc in the galactic disk. Two well-defined peaks 
of EW(CaT) at 19.1 kpc and 23.5 kpc are clearly evident. 
For the outer ring propagation velocity of 90 km it 
takes exactly 50 Myr to cross the distance between these 
two peaks in Fig. The angular difference between the 
peak positions of 6.6" should be easily detected. Now given 
the linear difference between the primary and secondary 
peak positions i?c a T is known from observations of a ring 
galaxy, one can calculate the propagation velocity of gas 
density wave by dividing i?caT on the value of tcaT- Or 
alternatively, if propagation velocity is known from in- 
dependent measurements or/and detailed modeling, one 
can calculate -RcaT, thus breaking new ground for the 
distance-independent estimates of a ring galaxy linear di- 
ameter. Hence, EW(CaT) may provide the means for a 
Hubble-constant-independent estimate of a distance to a 
ring galaxy. 

This method of evaluating the propagation velocity of 
gas density waves is most effective in the solar metallicity 
environments. For sub-solar metallicities however the pri- 
mary and secondary peaks in EW(CaT) become much less 
pronounced as shown by the open squares in Fig. 1101 for 
the Zq/2.5 metallicity. In fact, the primary peak almost 
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Fig. 10. The model radial profile of EW(CaT), as would 
be expected if the Cartwheel's disk had the solar metallic- 
ity (filled triangles) and z Q /2.5 metallicity (open squares). 

disappear for metallicities below z©/5. Furthermore, the 
absolute values of EW(CaT) generally decrease for pro- 
gressively lower metallicities (Mayya 1997). These diffi- 
culties make the observational detection of the primary 
and secondary peaks problematic in the sub-solar metal- 
licity environments. For the Cartwheel galaxy, the solar 
metallicity is expected only in the central parts. The outer 
parts are most probably of sub-solar (zq/3.5) metallicity 
and this method will not work. Unfortunately, heavy ele- 
ment abundance measurements of Bransford et al. |1998) 
indicate that individual star-forming knots in the rings of 
all ring galaxies in their sample have sub-solar metallici- 
ties in the range from zq/2.5 to Zq/5. We have to look for 
metal-rich ring galaxies in order to apply this method for 
determining the propagation velocity of gas density wave. 

5.2. Other spectral lines 

Some other prominent absorption lines might be used for 
estimating the propagation velocity of gas density waves 
in ring galaxies. Synthesized spectra of a single-burst stel- 
lar population (PEGASE2 or Starburst99) indicate that 
the equivalent widths of such optical absorption lines as 
Mg b, Na D, and some of strong Fe lines have a local max- 
imum/minimum between 10 and 250 Myr. A set of these 
lines can help to establish " time labels" along the radius of 
a ring galaxy under an assumption that most stars in the 
galaxy were born in a single star formation event. Even 
two localizations of such extremes in the radial profiles of 
absorption line equivalent widths can provide us with the 
ages of stellar populations at different radii and hence give 
us an insight into the value of propagation velocity of gas 
density wave that triggered the star formation event. Use 



of a large aperture telescope in combination with a con- 
temporary spectrograph makes it possible to obtain high 
signal-to-noise spectra of ring galaxies in a few hours of 
integration time (Moiseev[2002). 

5.3. Non-thermal radio emission 

Evolved stellar populations in the wake of an expanding 
ring density wave explode as supernovae. Type II super- 
nova remnants are responsible for strong non-thermal ra- 
dio emission. Condon and Yin (1990) argue that most 
Type II supernova remnants and their consequent non- 
thermal radio emission result from stars with masses just 
over 8 Mg, while stars more massive than 20 Mq ac- 
count for most of the ionizing radiation and consequent 
Ha emission. If so, we might expect to observe the dif- 
ference between the peak positions of the radial profiles 
of non-thermal radio emission and Ha surface brightness 
in ring galaxies. Indeed, the ring galaxy VII Zw 466 does 
show signs of such a difference when radio emission peaks 
on the inside edge of the disk and does not coincide with 
the prominent HII region complexes. The radial difference 
between the peaks in Ha and non-thermal radio emission 
might be used to infer the propagation velocity of gas 
density waves in ring galaxies in the same manner as in 
Sect. El 

6. Conclusions 

Our numerical hydrodynamics modeling of ring galaxy 
formation demonstrates that the maximum radial expan- 
sion velocity of gas (w gas ) in the first ring is always be- 
low the propagation velocity of the first ring itself (i> r ing) 
The inequality w gas < Wring is less strict in ring galaxies 
of comparable disk/halo masses than in halo-dominated 
ring galaxies. Numerical hydrodynamics modeling of the 
Cartwheel galaxy indicates the outer ring is currently 
propagating at u r i ng « 100 km s _1 , while the maximum 
radial expansion velocity of gas in the outer ring is cur- 
rently u gas Ri 65 km s _1 . The latter value is in marginal 
agreement with the measurements of Higdon (1996) based 
on HI kinematics, but substantially disagree with the mea- 
surements of Amram et al. 1)19981) based on kinematics of 
ionized gas. We conclude that the dynamical properties 
of neutral and ionized gas may differ significantly in ring 
galaxies. 

We demonstrate that the radial B — V/V — K color 
gradients, together with the Qbvk extinction-free indices, 
can be used for estimating the value of Wring m ring galax- 
ies. Modeling of the radial B — V/V — K color gradi- 
ents measured in the Cartwheel disk by Marcum et al. 
(1992) is sensitive to the assumed value of a Hubble con- 
stant. Nevertheless, the color gradient modeling imposes 
the lower limit of Wring — 40 km s _1 on the propagation 
velocity of the outer ring for any value of a Hubble con- 
stant in the range H — [50 — 100] km s _1 Mpc -1 . For the 
adopted in this work H Q = 65 km s" 1 Mpc" 1 , the best 
agreement between the model and observed radial color 
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gradients is achieved for velocities of 90 km s -1 and even 
higher. 

Our modeling of the Cartwheel galaxy demonstrates 
that the angular differences between the peaks in the 
Ha, K-, and 5-band radial surface brightness profiles are 
Ha - #=1.65" and Ha - 5=1.05", which should be de- 
tected with the achieved resolution. On the other hand, 
we find that the azimuthally averaged radial Ha surface 
brightness profile of the Cartwheel's outer ring does not 
peak exterior to those in Ks- and B-bands, contrary to the 
predictions of our model for v r ing = 90 km s _1 . We con- 
clude that a 41° inclination, along with a finite thickness 
and warping of the Cartwheel's disk, act to smear out the 
angular differences between the peaks in the Ha, Ks, and 
B radial profiles. Indeed, we construct a series of the ra- 
dial surface brightness profiles azimuthally averaged over 
a 1° sector. These pie-chart-like radial profiles do show 
an expected tendency and the Ha radial profile peaks ex- 
terior to those in Ks- and B-bands at the position of the 
Cartwheel's outer ring. The Ha — Ks and Ha — B peak dif- 
ferences are the biggest along the Cartwheel's major axis, 
where effects of inclination and finite thickness are mini- 
mized, and are in close agreement with our model predic- 
tions. The peak differences obtained along the Cartwheel's 
major axis imply Wring ~ HO km s . Azimuthal vari- 
ations in Ha — Ks, Ha — B, and B — Ks peak differ- 
ences support our assumptions about a thick and non- 
flat nature of the Cartwheel's outer ring. We conclude 
that the Cartwheel's outer ring is currently propagating 
at w r ing ~ 100 km s , which is in agreement with our 
numerical hydrodynamics results. 

We show that the equivalent widths of some spectral 
lines, such as the Call triplet, might be used for determin- 
ing the propagation velocity of ring density waves in so- 
lar metallicity environments. Solar metallicity rings have 
not been detected so far (Bransford et al. 1998), but they 
may be expected in young ring galaxies where the rings 
have not yet expanded out to the metal-unpolluted re- 
gions. Non-thermal radio emission is also a promising tool 
for determining the wave velocity in ring galaxies with 
high rates of star formation like the Cartwheel. 
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